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ABSTRACT 

We determine both representations of the Fundamental Plane [R c oc <Tq (/)~ 6 andR oc (<Tq (I) ( 7 1 ) a ] 
and the luminosity-effective phase space density (L oc / e ~ 7 ) scaling relation for A^-body remnants of 
binary mergers of spiral-like galaxies. The main set of merger simulations involves a mass-ratio of 
the progenitors in the range of about 1:1 to 1:5, harboring or not a bulge- like component, and are 
constructed using a cosmological motivated model. Equal-mass mergers are also considered. Remnants 
lead to average values for the scaling indices of (a) k 1.6, (b) rs 0.6, (A) rj 0.7, and (7} rj 0.65. 
These values are consistent with those of JsT-band observations (Mobasher et al. 1999) of ellipticals: 
(a) R5 1.5, (b) « 0.8, (A) rss 0.7, and (7) « 0.60. The b index is, however, not well reproduced. 
This study does not allow us to establish a conclusive preference for models with or without a bulge 
as progenitors. Our results indicate that the L-f c and FP scalings might be determined to a large 
extent by dissipationlcss processes, a result that appears to be in contradiction to other dissipationlcss 
results. 

Subject headings: galaxies: ellipticals - galaxies: kinematics and dynamics - methods: numerical, 
A-body simulations 



1. INTRODUCTION 

A large set of observational (e.g. Schweizer 1998, 
Struck 2005, Rothberg & Joseph 2004, 2006) and theoret- 
ical (e.g. Barnes 1998, Burkert & Naab 2003, Naab, Jes- 
seit & Burkert 2006) evidence shows that elliptical galax- 
ies are consistent with the picture of being formed by 
mergers of spiral galaxies (Toomre 1977). However, some 
matters are still open to discussion, such as the problem 
of stellar ages and metallicities (e.g. Peebles 2002, Ren- 
zini 2006, Naab & Ostriker 2007). 

Observationally, one global scaling relation displayed 
by ellipticals is the Fundamental Plane (FP; Djorgovski 
& Davis 1987, Dressier et al. 1987), a linear-log rela- 
tion among the effective radius R c , the central stellar 
1-dimensional velocity dispersion <7q, and the mean ef- 
fective surface brightness (2%: R c oc a^(I)^ b . Values 
of a and b depend on several factors among the most 
important are: (i) the observed wavelength, (ii) the fit- 
ting method, (hi) the assumed underlying light profile 
in early-type galaxies and, (iv) the width and brightness 
of the magnitude range of the sample [e.g. J0rgensen, 
Franz & Kjaergaard 1996, Pahre et al. 1998, Mobasher 
et. al 1999 (MGZA99), Bernardi et. al. 2003 (B03), Jun 
& Im 2008, D'Onofrio et al. 2008, La Barbera et al. 
2008, Nigoche-Netro et al. 2009]. A wide range of val- 
ues a Ri (1-2) and b ss (0.5-1) can be found in the diverse 
cited works. 

If structural and kinematical homology is assumed, as 
well as a constant mass-to-light ratio, direct application 
of the physical virial theorem to observational related 
quantities leads to a FP with coefficients (a, b) = (2, 1). 
The discrepancy of these values with the observed ones 
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is usually termed as the "tilt" of the FP; see also (|4]). 

The overall features of the observed FP are repro- 
duced by merger studies [e.g. Gonzalez-Garcia & van 
Albada 2003, Aceves & Velazquez 2005 (AV05), Boylan- 
Kolchin et al. 2005, Robertson et al. 2006 (R06), Dekel & 
Cox 2006] , but all show differences in values obtained for 
the indices a and b depending on the details of the initial 
conditions (e.g. the presence or not of a bulge and/or a 
gas component) used in the simulations. It is likely that a 
combination of broken structural and dynamical homol- 
ogy and stellar population effects, along with a variation 
of the central mass-to-light ratio, may lead to the expla- 
nation of the tilt (e.g. Scodeggio et al. 1998, Trujillo, 
Burkert & Bell 2004, AV05, Riciputi et al. 2005, R06, 
Bolton et al. 2007, Hyunsung & Im 2008). 

The study of La Barbera et al. (2008) points in the di- 
rection that the tilt of the FP is not driven by stellar pop- 
ulations but from other processes like a homology break- 
ing. They derive tight values of a w 1.5 and b w 0.75 
from a large sample of early-type galaxies by combining 
SDSS and UKIDSS data in the optical and NIR. How- 
ever, D'Onofrio et al. (2008) find larger variances of the 
distributions of the FP coefficients as derived from a sam- 
ple of three larger surveys (WINGS, NFPS and SDSS) 
suggesting that a such wide range of values for a and b 
seems to be in contradiction with the idea of a universal 
FP relation and speculate that the FP corresponds to a 
bent surface. 

The importance of dissipative processes in the merger 
hypothesis was early recognized from estimates of the 
central physical phase-space density of ellipticals, that 
turned out to be higher than that of spirals (e.g. Os- 
triker 1980, Carlberg 1986, Gunn 1987, Lake 1989, Ko- 
rmendy 1989). In a dissipationless scenario of galaxy 
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merging, Liouville's theorem (e.g. Binney & Tremaine 
1987) demands that the central physical phase-space 
density remains constant so, given the previous obser- 
vational estimates, it prohibits a dissipationless picture 
of merging disk galaxies to form ellipticals. Robert- 
son et al. (2006) using a large set of simulations of equal- 
mass mergers without and with a gas component, and 
other physical processes, conclude that the FP in the 
K band (Pahre et al. 1998) and the effective radius - 
stellar mass scaling relation [_R C -M»; Shen et al. 2003 
(S03)] can be only reproduced by mergers of disk galax- 
ies with a gas fraction / gas > 30% with respect to the 
mass of the disk component, M<j; ruling out dissipation- 
less mergers as mechanism to produce these two global 
scaling relations. Cox et al. (2006) also find that / gas is 
of the same order to reproduce the kinematic properties 
of ellipticals. 

Obtaining higher central phase densities in models of 
elliptical galaxy formation can be achieved, aside of tak- 
ing into account dissipational processes, by including a 
bulge-like component in the progenitors (e.g. Vedel & 
Sommer-Larsen 1990, Hcrnquist, Spergel & Heyl 1993; 
hereafter HSH). This result is in agreement with several 
TV-body simulations that have found that mergers of pure 
disk galaxies, residing inside dark halos with a core-like 
profile, are not able to reproduce the surface brightness 
density profiles typically found in ellipticals (e.g. Garcfa- 
Gonzalez & Balcells 2005, Naab & Trujillo 2006). 

The combination of arguments based on the non- 
reproducibility of the FP, the estimates of phase-space 
densities and the density profiles obtained from diverse 
TV-body simulations have lead to the idea that collision- 
less models of disk galaxies, with or without a bulge, are 
ruled out as viable progenitors to form ellipticals (e.g. 
Carlberg 1986, HSH, Mao & Mo 1998, R06, Ohorbe et 
al. 2006, Robertson et al. 2006). Nonetheless, there are 
results of iV-body mergers of pure disk galaxies residing 
in cuspy dark halos that lead to density profiles consis- 
tent with those of early-type galaxies (Aceves, Velazquez 
& Cruz 2006: AVC06) and are able to reproduce ade- 
quately the Fundamental Plane scaling relation (AV05). 
Our results of dissipationless mergers (see below) suggest 
that they can not be ruled out on the basis of scaling ar- 
guments. 

On the relevance of a bulge-like component in progen- 
itors of merger remnants, there is observational evidence 
that an important fraction of spirals do not harbor a 
bulge-like structure, or that it does not have a signif- 
icant contribution, as is the case of late type morpho- 
logical Hubble types usually called "flat galaxies" (e.g. 
Goad & Roberts 1981, Karachentsev 1989, Kautsch et al. 
2006, Dommgucz-Palmero et al. 2008). According to the 
sample considered by Kautsch et al. about one third of 
spirals are flat galaxies or simple disks, and in the sample 
of Dommguez-Palmero et al. about 30% (54) harbor a 
bulge and about 70% (137) have no measurable bulge. 
Thus, these flat galaxies arise the immediate issue about 
to what extent A-body simulations are able to establish 
global scaling relations like the FP when merging. 

Undoubtedly, an important test on the details of the 
merger hypothesis is imposed by the phase-space con- 
straints. A comparison among theoretical models of the 
central physical phase-space density, / p , using different 
models and initial conditions, although important for our 



understanding of the structure of remnants, needs to be 
related to observations in order to determine their rele- 
vance. Different works have already shown that f p in- 
creases with the inclusion of a dissipative process (e.g. 
R06) or a bulge inside a core-like halo (e.g. Naab & 
Trujillo 2006). However, given that the physical central 
phase-space density, / p , cannot be determined observa- 
tionally an estimator has to be used when comparing 
with observations. 

An early estimator was the core coarse-grained phase- 
space density, / c , proposed by Carlberg (1986). How- 
ever, HSH considered that a more robust one is the ef- 
fective coarse-grained phase-space density, f c . A strong 
dependency of f c with the luminosity, L, of ellipticals 
has been found in several works (e.g. Carlberg 1986, 
Lake 1989, HSH, Mao & Mo 1998). For example, using 
data of Bender, Burstein & Faber (1992), HSH found 
that Lb oc / e -0 ' 54 for the compact, intermediate and gi- 
ant ellipticals with a RMS scatter of 0.15 mag; with Lb 
being the total luminosity in the £?-band. 

The FP and L-f e scaling relations appear to be tight 
enough to serve as constraints for the formation scenario 
of elliptical galaxies. However, it should be said that up 
to now there is no direct comparison of the L-f c scaling 
relation obtained from observations with its counterpart 
provided by the remnants of merger simulations. Pre- 
vious works (e.g. HSH, R06) have computed only the 
physical cumulative coarse-grained distribution function, 
s(f), which is compared with a model of an elliptical 
galaxy or cither with the effect of introducing a bulge 
and/or gas component in the progenitors, but it has not 
been related with an observational quantity. 

By restricting the comparison of our numerical rem- 
nants with infrared observations, for which stellar popu- 
lation age effects are minimum (e.g. Bruzual & Chariot 
2003), and by adopting a procedure similar to the one 
used in observational studies, we examine to what ex- 
tent A-body numerical simulations of spiral-like galaxy 
mergers, with or without a bulge-like component, are 
able to reproduce the observed L-f e , the FP and the R- 
global relations. The latter two results are compared 
with those obtained in our earlier works (AV05, AVC06) 
where the progenitor galaxies lacked a bulge component 
and with recent results from other authors. 

It is important to mention that we do not attempt 
to determine the "zero-point" of the above scaling rela- 
tions since this surely requires to include more physical 
processes, such as star formation, cooling, and feedback 
among others; a matter that is out of the scope of the 
present work. 

This paper has been organized as follows. In Section 2 
we summarize the numerical galaxy models and tools 
used to carry out this work. Section 3 presents our results 
and Section 4 contains a discussion and final comments. 

2. SIMULATIONS 

The method to build up our galaxy models is described 
in AV05 and for the sake of clarity it is briefly sum- 
marized. The galaxy models considered here satisfy the 
Tully-Fisher relation (Tully & Fisher 1977). They are set 
up by combining the Press-Schechter formalism of hierar- 
chical clustering and the cosmologically motivated model 
of disk galaxy formation of Mo, Mao & White (MMW; 
1998). Essentially the method outlined by Shen, Mo & 
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TABLE 1 

Properties of bulge progenitors of mergers. 
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Mrf, Rd, Zd and correspond to the mass, radial scale length, vertical scale length and number of particles for the disk component. M/j, 
7"200; -V c an d indicate the mass, virial radius, spin parameter and concentration of the progenitor halos. The last two columns refers 
to Eulcr's random angles (in rads) used to rotate the axi-symmetric galaxy. In all cases, the bulge component, represented by a Hernquist 
profile, was constrained to have a mass, scale length and number of particles of M^/4, Rd/3 and N^/i, respectively. 



Shu (2002) is followed. In this scenario a formation red- 
shift for disks needs to be established, which was set here 
to be z — 1 corresponding to a look-back time of w 8 Gyr 
(Peebles 1993); simulations with disks formed at higher 
z were not considered. In the MMW model, disks at 
higher redshifts are smaller and denser than those at the 
present cosmic epoch. A ACDM cosmological model with 
TO = 0.3, £!a = 0.7 and a Hubble parameter h — 0.7 is 
adopted. The numerical progenitors satisfy the disk sta- 
bility criteria of Efsthatiou, Lake & Negroponte (1982), 
v^ n /(GR^ 1 Md) > 1.0 where v m , Rd and Md are the max- 
imum rotational velocity, the radial scale length and the 
total mass of the disk. 

The iV-body galaxy models consist of a spherical dark 
halo, a stellar disk and a bulge component. The dark halo 
follows a Navarro, Frenk & White (1997) profile with an 
exponential cutoff. The disk component has an expo- 
nentially decaying radial density profile and an isother- 
mal vertical structure. Finally, the bulge-like component 
is represented by a spherical Hernquist profile (Hern- 
quist 1990). The bulge in this galaxy formation scenario 
is introduced by following the procedure described by 
Springel & White (1999). In all fixed mass frac- 

tion of 0.25Md and a scale length of one-third of the 
radial scale of the disk are adopted for the bulge. Table 
1 summarizes the parameters of our progenitor galaxies 
with a central bulge component. 



Also, for comparison purposes the remnants of our pure 
disk galaxy models described in AV05 together with five 
new similar merger simulations are included. In this case, 
the progenitors of our remnants resemble the " flat galax- 
ies" studied by Kautsch et al. (2006) and Dommguez- 
Palmero et al. (2008). 

Only parabolic encounters are considered with pericen- 
ter parameters (5-20 kpc) consistent with those found in 
cosmological simulations (e.g. Navarro, Frenk & White 
1995). The spin orientation of each galaxy, relative to 
the orbital plane and defined here using Euler's angles 
(Goldstein 1950), has been taken randomly. We have 
only considered binary encounters, so our results may 
not apply to massive early-type galaxies since they are 
probably the result of multiple disk galaxy merger events 
and/or mergers of early- type galaxies (e.g. Weil & Hern- 
quist 1996; Naab, Khochfar & Burkert 2006). 

The progenitor galaxies of our main set of merger sim- 
ulations have a mass-ratio in the range of about 1:1 to 
1:5. However, we have also added a new set of ten equal- 
mass mergers to assess the importance of the progen- 
itors masses in dissipationless mergers of disk galaxies 
in determining the FP and the L-f e global scaling rela- 
tions. These simulations were obtained using as progen- 
itors galaxies the ones partaking in mergers M01, M05, 
M07, M08 and M10 of AV05 that lead to luminous- 
dominated merger cores. New orbital parameters, fol- 
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TABLE 2 

Remnants properties with bulge progenitors. 
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Mbl2 


2.58 


8.55E+00 


-2.12E+01 


-2.78E+01 


1.54E+02 


7.90E+00 


1.6IE+02 


-2.70E+01 






1.10E+01 


-2.05E+01 


-2.77E+01 


1.56E+02 








Mbl3 


2.86 


5.14E+00 


-2.18E+01 


-2.74E+01 


1.30E+02 


5.07E+00 


1.31E+02 


-2.66E+01 






6.37E+00 


-2.14E+01 


-2.74E+01 


1.32E+02 








Mbl4 


2.19 


3.04E+00 


-2.32E+01 


-2.76E+01 


1.98E+02 


3.68E+00 


1.78E+02 


-2.70E+01 






6.01E+00 


-2.21E+01 


-2.80E+01 


1.94E+02 








Mbl5 


4.02 


3.95E+00 


-2.24E+01 


-2.74E+01 


1.43E+02 


4.07E+00 


1.38E+02 


-2.66E+01 






3.98E+00 


-2.24E+01 


-2.73E+01 


1.43E+02 









First line of each merger corresponds to a Sersic fit, while the second to a i? 1 / 4 . Simbols n, R c , (I)c, Ml, and ctq correspond, respectively, 
to the Sersic index, effective radius, mean effective surface brightness, total "magnitudes" (—2.5 log Ml), and central velocity dispersion 
within R c /S. Values R50, ""50 an( i ^h,50 are the proyected half-luminous mass radius, the ID velocity dispersion, and the half-mass 
proyected luminous mass (in "magnitudes"), respectively, of our remnants measured directly from the simulations. Units are indicated. 



lowing the above procedure, were obtained. 

Simulations were done using the parallel tree-based 
code Gadget-2 (Springel, Yoshida & White 2001, 
Springcl 2005). They were evolved for about 8Gyr, a 
time similar to the time spanned from z = 1 to the 
present epoch and at which the remnants had reached 
equilibrium. Softening for disk, bulge and halo particles 
were taken to be td — e& = 35 pc and th — 350 pc, 
respectively, with a parameter ErrTolForccAcc= 0.0025. 
Energy conservation was better than 0.75% in all simu- 
lations. 

3. RESULTS 

In this section we present the FP, R c -M* and the L- 
/o relations for our dissipationlcss mergers of spiral-like 
galaxies. 

We resort to a procedure similar to the one used in 
practice with the observations in order to avoid, as far as 
possible, any kind of biases in our results difficulting their 
comparison with them, if-band observations are better 
suited for our purposes, since they are less affected by 
stellar population effects. For this reason we have cho- 
sen the observational data of MGAZ99 which is publicly 
available. 

The results for bulgeless progenitors are obtained by 
considering only those mergers of our sample (10 out of 
15) that resulted in remnants dominated by "luminous" 



matter inside their effective radius (R < R c m 1-2 kpc) 
as is suggested by observations (e.g. Gerhard et al. 2001, 
Thomas et al. 2007). However, the scaling indices do not 
change significantly if remnants centrally dominated by 
dark matter are included. 

Each remnant was "observed" along 100 random dif- 
ferent lines-of-sight and fitted with a R 1 ^ 4 (de Vau- 
coleurs 1953) and a Sersic profile (Sersic 1968). The 
was adopted to compare directly with 
MGAZ99 while the Sersic profile is better suited to fit 
our numerical remnants and it is characterized by the 
following surface density: 

S(i?) -EocxpRWRe) 1 /"], (1) 

where b = b(n) w 2n - 1/3 + 4/(405n) (e.g. Graham & 
Driver 2005), S is the central projected surface luminous 
mass density and n is the corresponding Sersic index. 
The radial interval of the fit is taken from our numerical 
resolution value for the luminous matter (2.8ed = 2.8efc w 
100 pc) to the outer radius enclosing 95% of the projected 
luminous mass, determined directly by counting particles 
from the simulations. The parameters R e and n of 
the profile were found by a ^-minimization using the 
Levenberg-Marquardt method (e.g. Press et al. 1992). 
The total luminous mass of a Sersic model is given by 

M L = ^r(2n)S i? c 2 , W 
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TABLE 3 
Scaling indices. 
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were T is the Gamma function. For a crude estimate 
of the luminosity, L, of the numerical remnants we may 
choose a constant stellar mass-to-light ratio, T*, leading 
to L = T^Ml. To estimate a total magnitude Mx in 
the if-band we take as a fiducial values T„ = 0.97, in 
solar units, and M Q = 3.32. 

The central projected, luminosity weighted, velocity 
dispersion, <jo, of the luminous particles was computed 
inside a circular region of radius i? e /8 for each remnant 
projection, as is usually done in observational studies 
(e.g. J0rgensen et al. 1996, MGAZ99, D'Onofrio et al. 
2008). A value for the mean effective surface brightness 
(7) e can be obtained from the fitted parameters of the 
profile and by assuming a constant T*. Results of dif- 
ferent quantities for our set of simulations with a bulge 
progenitors are shown in Table 2. 

3.1. Fundamental Plane 

A log-plane of the form: 

log R c oca log cr + b log(I) e (3) 

was fitted to our remnants values by two standard meth- 
ods (e.g. B03): the direct least-square method (lsq) 1 
and an orthogonal plane (ort) procedure. The tilt A of 
the FP given by 

log J R o cxAlog( ( T 2 (/)- 1 ), (4) 
was also determined using the two previous fitting meth- 
ods. Both the LSQ and ORT methods were tested against 
synthetic data, in 2-D and 3-D, with an excellent agree- 
ment between the synthetic and fitted parameters. 

Table 3 summarizes the main results for our simula- 
tions with and without a bulge component, and using 

1 Formulae (9) and (14) in Robertson et al. (2006) for the di- 
rect least-square method display two misprints; compare them with 
Bernardi et al. (2003). 



either a Sersic (S) or i? 1 / 4 profile. Columns (3-5) pro- 
vide the indices a and b along with the corresponding 
RMS. Column (6) provides the tilt (A) of the FP, and 
Column (7) shows the index /i for the fitting of the re- 
lation R oc M*. Column (8) presents the index 7 of the 
luminosity- phase space relation (see below), and the last 
column (9) indicates the fitting method used. We have 
included also results for a set of simulations of equal-mass 
mergers (eM) lacking a bulge component. Errors in the 
ORT fits were estimated by bootstrap, while those in the 
LSQ fits is the standard RMS. 

For comparison purposes, we have included in Table 3 
the following: (a) values quoted by Robertson et al. 
(2006) for their N-body and full (iV-body+gas) simu- 
lations. The methodology followed by these authors is, 
however, different to the one adopted here and in ob- 
servational studies. Robertson et al. measure directly a 
half-mass stellar effective radius (R50) inside which the 1- 
dimensional dispersion velocity, 050, and a mean surface 
density I e — Ml,5o(< R^o)/ttR\q are determined; a 3D 
log-plane ([3]) is then fitted by a LSQ method, (b) Some 
fits from observational data are also included. Data from 
MGAZ99 for Coma ellipticals are used to compute differ- 
ent scalings. We indicate Pahre et al. (1998; PDdC98) 
average fitted parameters for the FP for all the clusters 
in their sample as well as the corresponding ones to the 
Coma cluster. The FP in the K-band from La Barbera et 
al. (2008; LaB08) is also shown. Despite the SDSS data 
are in the optical region, we include the FP and i?-M* 
relations derived by Bernardi et al. (2003) and Shen et 
al. (2003) using a i? 1 / 4 -LSQ-ORT and a Sersic-ORT profile, 
respectively. 

In Figure [1] we plot the FP in the "virial plane" repre- 
sentation: Re-a^I)^ 1 . Here, results from our iV-body 
mergers, with (solid dots) and without (filled triangles) 
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Fig. 1. — Comparison of near-infrared (NIR) Fundamental Plane 
(Mobasher et al. 1999, open squares) with our TV-body remnants, 
with (filled circles) and without (filled triangles) a bulge. Results 
for equal-mass mergers, without a bulge (NB eM, open stars), are 
also shown. The "tilt" A of the FP in each case is indicated: (a) 
mergers with (long-dashed line) and without (solid-line) a bulge 
component, (b) equal-mass mergers without a bulge (dotted line) 
and, (c) K-band observation by MGAZ99 (dot-dashed line). 

a bulge-like component, are shown. Equal-mass mergers 
lacking a bulge are also included (open stars). For com- 
parison, the FP in the K-band is represented by open 
squares and, finally, the tilt, A, in each case is provided. 

From Table 3, it follows that our iV-body remnants 
show different values for the indices a, b, A and /x, de- 
pending on both the fitting method and the assumed 
profile for the luminous matter. Averaging over all the re- 
sults, ignoring the details of the fitting, lead to (a) fts 1.6, 
(b) w -0.6, (A) w 0.7, and (/x) « 0.5. On the other 
hand, the corresponding K-b&nd averages are: (a) ps 1.5, 
(b) re -0.8, (A) « 0.7, and (fi) « 0.6. 

The previous average values indicate a general good 
agreement between the theoretical and observational val- 
ues; they are also consistent with those obtained by 
Bolton et al. (2007) using lens galaxies: a = 1.50 ± 0.32 
and b = —0.78 ± 0.13. However the index b, associated 
with the luminosity and not the stellar mass, is the one 
that shows a larger discrepancy from the observed val- 
ues. A similar thing occurs for the fi index. Nonetheless 
the tilt, A, obtained through the fit of the combined vari- 
able <7o/(/) c , yields results that agree very well with the 
observational data. 

Judging only from the individual indices a and b it is 
not clear if models of progenitors with or without a pri- 
mordial bulge are to be preferred, although progenitors 
without a bulge and different mass-ratios appear to be 
favored when the tilt A is estimated. However the rather 
small statistics of our study precludes a conclusive an- 
swer. 

3.2. Luminosity - Phase Space Relation 

The effective coarse-grained phase-space density of our 
numerical remnants is computed using equation (2.12) of 



HSH for the luminous matter, namely: 

/c = -W > ( 5 ) 

where the gravitational constant G has been taken to be 
unity. This estimator allows a more direct comparison 
of the global trend of the coarse-grained central phase- 
space with Mk of our numerical remnants with their 
observational counterparts. A direct least-squares fit of 
the form 

\ogL K oc -7 log/ c , (6) 

where 7 is the fitting parameter, is done to the mean of 
the values of the projections of each remnant. The mag- 
nitude value is taken as Mk — —2.5 logi^- A similar fit 
is done to the observational data of MGAZ99. 

The correlations M^-fe for our numerical remnants 
corresponding to progenitors with (solid dots) and with- 
out a bulge component (solid triangles) are indicated in 
Figure [5J Results for the equal-mass mergers are shown 
with an open-star symbol. For comparison, the Mx—fe 
relation derived from the data of MGAZ99 is also shown 
(open squares). The index 7 obtained for the different 
cases considered in this work are shown in Table 3. 

A scaling relation of the form L cx //T 7 results from 
these fits, with 7 = 0.66 and 7 = 0.51 for our remnants 
(using a i? 1 / 4 profile and lsq) from progenitors with and 
without a bulge, respectively. For equal-mass mergers, 
we have 7 = 0.63. 

The value for the MGZA99 sample is 7 = 0.57. This 
value is similar to that obtained in the B-band for the 
data of Bender et al. (1992) by HSH (7 = 0.54), and that 
if we use the ellipticals from D'Onofrio (2001) (7 = 0.51); 
despite the fact that D'Onofrio uses a Sersic profile to fit 
the brightness distribution while Bender et al. consider 
a i? 1 / 4 profile. 

In general good agreement, within uncertainties, 
among the iV-body and observed 7 indices exist. Taken 
at face value, non-equal mass (bulgeless) mergers have a 
closer value to the observed ones. However, as for the FP 
case, our results preclude us from favoring a particular 
model for the progenitors. 

4. DISCUSSION AND FINAL COMMENTS 

An important conclusion of recent works on mergers 
(e.g. R06, Naab & Trujillo 2006, Cox et al. 2006, Dekel 
& Cox 2006) is that dissipation plays a major role in 
reproducing different scaling relations; such as the FP 
or the R e —M* relation. In particular, R06 find that the 
initial fraction of gas in merging galaxies should be ~ 
30% of the disk mass in order to reproduce both the 
observed scaling of the Fundamental Plane and the R c - 
M* relation. 

However, previous works by AV05 and AVC06 were 
able to reproduce rather adequately these global scaling 
relations without invoking a gas component or even a 
bulge-like structure. The b index of the FP was less well 
reproduced, perhaps due to the assumption of a constant 
Y* to obtain (J) e . The results found here for progenitors 
with a bulge, and bulge-less equal-mass mergers, are in 
general in same standing when compared with NIR data 
under methods that mimic the observational procedure. 

As mentioned earlier, the observed tilt A is well re- 
produced in our dissipationless simulations. The reason 
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Fig. 2. — Effective coarse-grained phase-space density / e versus 
total magnitude My. The it-band data sample of ellipticals from 
MGAZ99 (open- squares, dot-dashed line) yields a scaling index of 
7 = 0.57. The scalings for our remnants with a bulge (long-dashed, 
solid circles) and without it (solid line and triangles) are shown; 
the bulge- less equal mass merger are also plotted (dotted line and 
open stars). 

for this agreement, although the b index is rather poorly 
matched (see Table 3), might be related to correlations 
present among the variables (and their errors) involved 
in the fitting of the FP. Such correlations will affect the 
fitting parameters depending on how the FP is fitted, 
either as in equation © or in (|3]). 

On the other hand, it is not clear the origin of the dif- 
ference between our fitted values of a and b of the FP and 
those derived by Robertson et al. (2006) for their equal- 
mass dissipationless simulations. For instance, when we 
use the half-projected stellar mass radius, R50, and the 1- 
D velocity dispersion, 050, inside this radius we obtain a 
value of a w 1.8 that tends to the expected virial value as 
in R06. However, even in this case our value for b sa 0.5 
is far from unity (see Table 3). The tilt A « 0.8 in this 
case tends to increase a little but still it is not unity. 

Given our results, the FP appears to arise mainly from 
broken structural and dynamical homology between the 
observed (luminous component) quantities {ao,R Cl (-Oe} 
and their theoretical counterparts {V 2 , Rq, M}; where 
V 2 = 2T/M, R G = GM 2 /\U\, and M is the total bound 
mass of the system. Here T is the total kinetic energy of 
the system, with contributions from rotational and ther- 
mal motion in remnants. An indication of the previous 
was presented in Aceves & Velazquez (2005b) . Variations 
in the stellar mass-to-light ratio, Y*, can be responsible 
for some of the tilt in the FP, especially if we focus on 
the value of the b index for our remnants. However, the 
effect of T* would tend, according to our results, to di- 
minish when the A parameter is used as an indicator of 
the FP tilt. 

The L-f e scalings shown in Figure 2 are in the same 
direction as the preceding results on the FP, suggesting 
that the establishment of this scaling relation is consis- 
tent with a dissipationless merging scenario. It is likely 



that dissipative processes (e.g. star formation and feed- 
back) will play a major role in defining the zero-point of 
all these scaling relations, and affect the indices in a not 
simple way. 

From energy arguments in a dissipationless scenario, 
and assuming equal mass galaxies, Hernquist et al. (1993) 
derived the particular relation L cx / c -0 ' 5 . When non- 
identical progenitors are taken into account these energy 
arguments predict a shallower log-slope; for example, for 
a mass ratio of about 1/5 it leads to L oc / c ~ ' 37 - So, 
according to these energy arguments we should expect a 
range of values in the interval 0.37 < 7 < 0.5 for the non- 
equal mass ratios considered here. Thus, the theoretical 
expectation on the basis of these energy arguments is 
about 30% lower than the average value ((7) « 0.7) for 
our no-bulge equal-mass mergers. Meanwhile these theo- 
retical values are more closer to the ones of our remnants 
from bulgeless progenitors. 

The average result from all of our mergers yield an 
index (7) w 0.65 that is in good agreement with the 
observed NIR value 7 w 0.60. This agreement spans 
for about four magnitudes and compromises remnants 
that can be cataloged from bright dwarf ellipticals to 
intermediate early-type galaxies (see Figure 2). 

It should be noticed that differences arise when a pri- 
mordial bulge is considered or not, and merger remnants 
from bulgeless progenitors tend to produce a scaling more 
close to the NIR one. In general, for larger "luminosity" 
the effective phase-space density tends to a higher value 
for remnants with progenitors having a bulge as is found 
in other works that quantify the physical distribution 
function f p (e.g. HSH, Naab & Trujillo 2006, Robert- 
son et al. 2006). At low "luminosity" the behavior of 
/ e does not show a clear difference for bulge or bulgeless 
progenitors. 

A physical explanation for the details of the behavior 
of L—f c (or it physical counterpart) is non-existent at 
the moment, although the arguments of HSH yield light 
on the subject. On other hand, unfortunately, the more 
extended study carried out by R06 did not address the 
L—f c scaling relation to compare directly with our N- 
body results. 

In summary, the observed A-tilt of the FP and the 7- 
index of the luminosity-phase space relation, both in the 
NIR, are consistent with those obtained here through 
dissipationless merger simulations; at least over the lu- 
minous mass range covered by our remnants. We ascribe 
such concordance to the use of progenitors constructed 
under a cosmological motivated model of galaxy forma- 
tion, with a NFW dark halo and satisfying initially a 
Tully-Fisher relation. 

As a concluding remark, given the results of this work, 
we would like to point out that trying to deduce the 
amount of primordial gas (/ gas ) in models of progenitors, 
taking part in merger events, to satisfy a global scaling 
relation like the FP, might be subject to a non-negligible 
uncertainty. It is likely that the zero-point of the scaling 
relations would be a better constraint to the dissipational 
component of disk galaxy progenitors. In our opinion this 
is a subject that requires further work, but it is out of 
the scope of the present work. 
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